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SUMMARY 

— ' The  two  essential  ingredients  oi  any  boundary  value  problem  ere  the  field 
equations  which  describe  the  physics  of  the  problem  and  a ret  of  relations  vhich 
specify  the  geometry  of  the  problem  domain.  Mesh  generators  or  grid  generators 
are  preprocessors  which  decompose  the  problem  domain  into  e lerge  number  of 
interconnected  finite  elements  or  curvilinear  finite  difference  stencils.  A num- 
ber of  such  techniques  have  been  developed  over  the.  past  decade  to  alleviate  the 
frustration  and  reduce  the  time  involved  in  ..he  tedious  manual  subdividing  of  e 
complex-shaped  region  or  3-D  structure  into  finite  elements «,  Our  purpose  hare 
is  to  describe  how  the  techniques  of  bivariate  and  trivariate  ^"blending  function"' 
interpolation,  which  were  originally  developed  for  and  applied  to  geometric  pro- 
blem* of  computer-aided  design  of  sculptured  surfaces  and  3-D  solids,  can  be 
adapted  and  applied  to  the  geometric  problems  of  grid  generation.  In  contrast 
to  ether  techniques  which  require  the  numerical  solution  of  complex  partial  dif- 
ferential equations  (end,  hence,  a greet  deal  of  computing),  the  transfioito 
methods  proposed  herein  are  computationally  inexpensive. 

1.  INTRODUCTION 

Over  the  past  decade,  a number  of  scheme*  have  been  developed  fer  automating 
the  generation  of  finite  element  and  curvilinear  finite  difference  grids.  Among 
these,  the  tttnsflai te  supping  technique  of  Gordon  and  Halil 3]  has  been  shown  to 
have  a number  of  advantages  icf.  !6},(7]).  Some  of  those  are: 

1.  Exact  modeling  of  boundaries 

2.  Minimal  input  effort 

3.  Automatic  node  connectivity  definition 

4.  Well-suited  to  interactive  graphic*  implementation 

5.  Good  correlation  botwatn  boundary  nodes  and  intar lor  mesh 

6.  Computationally  efficient 

7.  Esey  extension  to  three  dimension*. 

Wg  use  the  tens  "trans  finite"  to  do ter lb*  this  class  of  techniqusc  since, 
unlike  classical  methods  of  higher  dimensional  interpolation  which  match  the 
primitive  function  r it  i finite  number  of  point* , the  t ran* finite  methods  match 
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F at  a nondanumerable  number  of  points.  In  fact,  as  we  shall  see  below,  trane- 
finite  mappings  of  the  plane  match  F along  entire  curve  segments,  while  trans- 
finite  mappings  in  Euclidean  3-space  can  match  F exactly  on  the  six  faces  of  a 
curvilinear  parallelpiped. 

To  begin,  we  recall  the  geometric  interpretation  of  the  graph  of  a . ector- 
valued  function  of  two  independent  variables  s and  t 

F(e,t:  - [x.  U,t)  ,x,(s,t)  ,...,x  (s,t)]T.  (1) 

l i n 

2 

As  the  variallec  s and  t range  over  a dorn/.n  S in  the  3,t-plane  R , F(s,t) 
tra-ea  out  e region  R in  Euclidean  n-cpace.  En  That  is,  F maps  regions  in 
,R2  into  regions  in  En 

Fs  k2  -*■  En.  (2) 

For  two-dimensional  problems,  we  shall  be  concerned  with  continuous  transforma- 
tions F which  map  the  unit  square  S » ( 0,11x1 0,1]  one-to-one  onto  a simply  con- 
nected, bounded  region  R in  E2  or  E3.  Such  mapa  can  be  thought  of  as  topologi- 
cal distortions  of  the  planar  region  S onto  tho  two-dir'ensional  manifold  R, 

2 

which  is  ci tiier  a planar  region  (R  r-  S ) or  i surface  eabeddoO  in  3-epace 

(R  c E3}.  xn  either  case,  a ono-to-one  'un, valent)  wapping  S ♦ R ie  equivalent 

to  the  introduction  of  a curvilinear  co-ordinate  system  on  R.  The  curve  of  con- 

♦ _ 
etant  yenaralited  co-ordinata  a « »•  is  the  image  F(e*,t)  of  the  line  s»a*inS. 

Similerly,  the  curve  r(s,t»)  is  the  eat  of  all  points  in  R with  generalized 
co-ordinate  t ■ t*.  Thue,  the  point  F(e*,t*)  on  R ie  said  to  have  generalized 
co-ordinatee  ie»,t»),  end,  eince  the  sapping  S ♦ R is  univalent,  any  point  R 
can  be  uniquely  referenced  by  its  generalized  co-ordinetee . 

If  S i*  the  unit  cube  [0,l)x(0,l }x( 0 , 1]  in  the  a,t,u-parametar  space  R3  and 
R is  a bounded  region  in  Euclidean  3-cpace,  than  a one-to-one  mapping  f of  S on- 
to R can  be  envisioned  ae  a topological  distortion  of  the  cube  Into  R.  Such  e 
mapping  of  8 3 ♦ E3  generate.)  a curvilinear  co-ordinatization  of  the  solid  R so 
that  each  point  of  R may  be  referenced  by  its  gmnarmlixmd  coordinates  (s,t,u). 

For  bounded,  eiaply  connected  planar  domains  R,  one  could  of  course  generate 
en  orthogonal  co-ordinatization  by  means  of  e conformal  tapping  of  R onto  a 
canonical  region  such  ae  a square  or  a circle.  .Vc waver,  from  a practical  stand- 
point, the  construction  of  a conformal  map  ia  equivalent  to  the  solution  of 
Laplaca'a  equation  and  ia  thus  contrary  to  the  qoal  of  computational  simplicity. 
In  contrast,  the  traneflnite  napping*  described  below  are  relatively  simple  to 
construct  end  implement  for  a wide  variety  of  regions,  end  are  computationally 
inaxpansive. 


173 


2 . TWO-DIMENSIONAL  REGIONS 

We  first  consider  the  case  in  which  S is  the  unit  square  [0,l]x[0,l],  Let  us 
postulate  the  existence  of  a primitive  function  F which  maps  S onto  R.  It  should 

4 

be  remarked  that  the  'unction  f is  a 'iction  which  we  introduce  only  'or  nota- 

tional  simplicity  and  convenience.  In  practice,  the  only  thing  we  are  given  is 

the  geometric  description  of  R in  terms  of  its  boundary,  and  it  is  the  task  of 

the  analyst  to  cast  this  information  into  a form  appropriate  to  the  mapping 

formulas  considered  below.  This,  however,  is  not  difficult  to  do  and  can  be 

implemented  by  a computer  subroutine. 

*22 

Generically,  Ft  R ♦ E should  be  thought  of  as  a continuous  vector-valued 
function  of  the  two  independent  variables  s and  t such  that  3S  ■+  3R. 

For  example,  consider  the  following  mapping: 


?(s,t) 


4 st (l-s) Cl-t)  + (l+t/ZT) 


(l+t/^2)  sin  y- 


(3) 


This  maps  the  unit  square  [0,l]x[0,l]  onto  the  quarter  annulus  R shown  in  Fig.l. 
The  perimeter  of  the  unit  square  maps  onto  the  perimeter  of  R,  and  lines  of  con- 
stant s and  constant  t map  onto  the  curvilinear  co-ordinate  system  illustrated. 
In  other  words,  each  of  the  curves  shown  in  the  figure  is  a curve  of  generalized 
co-ordln*te  s • const,  or  t * const. 

-4 

Our  problem  is  to  construct  a univalent  (one-to-one)  function  U:  S ♦ R which 
watches  F on  the  boundary  of  S,  i-e. 

U(0,t)  - £(0,t),  U(s,0)  - F(s,0)  \ (4) 

3(1, t)  - F(l,t) , U(a,l>  - F(a,15  1 

a function  U which  interpolates  to  P at  a non-danuserable  set  of  points  as  in 
(4)  will  be  tensed  a trena  finite  interpolent  of  F. 

To  explain  the  xstion  of  transflnite  mapping,  vs  shall  find  it  convenient  to 
rely  on  the  algebraic  theory  of  approximation  developed  in  til  and  (4l.  in  this 
paper,  by  a projector  P we  shall  mean  an  idempotsnt  linear  operator  whose  do- 
main is  the  linear  space  F of  all  continuous  functions  da fined  on  S and  whose 

range  Is  a subspace  of  F.  The  above  interpolation  problem  (4)  can  be  viewed  as 

•*  * 

a search  for  a projector  P such  thet  U • P(r]  is  a univeisnt  sap  of  S * R which 

*4  *4 

satisfies  the  desired  interpoletoty  properties.  1)  is  termed  the  projection  of  F 
or  the  image  of  F under  P. 

Suppose  now  that  and  era  four  univariate  functions  which  satisfy 

tbs  cardinality  conditions 
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*l‘kl  ■ \k 5 

j 1/  i - k 

for  i,k  ■ 0,1 

to,  i + k 

' ) 

(S) 

V11  -*M 

for  1,1  ■ 0 , 1 ^ 

and  consider  the  projectors  P and  P defined  by 

fl  t 

P8[F)=40<8>F<s0,t)  + (^1(s)F(sl(t) 

Pt[F]  =^0(t)F(a,t0)  + ^(tJFCs,^) 

Then,  the  product  projection 
* 1 1 

* [ I ♦.(•)iK(t)FU  ,t.)  (7) 

8 C i-0j«0  1 3 X 3 

interpolates  to  F at  the  four  corners  of  [0,l]x[0,l]  and  the  Boolean  sum 
projection 


(Ps®  P )(F]  = P [FJ  + PjF]  - PsPt[F]  (8) 

interpolates  to  F on  the  entire  boundary  of  [0,l]x[0,l].  These  properties  of 
the  functions  (7)  and  (9)  may  be  readily  verified  by  evaluating  the  right-hand 
sides  for  the  appropriate  values  of  a and  t and  recalling  the  cardinality  pro- 
perties (5) i see  also  (3]  or  (4). 

The  functions  and  ^ in  the  above  formulae  are  as  yet  unspecified  except 
for  their  values  at  the  points  sQ  » tQ  •>  0 and  s^  » tj  ■ 1,  They  are  counonly 
referred  to  as  'blending  functions'  and  the  function  U ■ (P  • P )IfJ  is  termed 
a blended  interpolant.  The  simplest  choice  for  the  blending  functions  in  (5) 
is  the  set  of  four  linear  functions 

t 


v> 

$x(e> 


1 - s,  ^0<t)  - 1 
s.  ^(t)  - t 


(9) 


The  vector-velued  bivariate  function  0 3 (P  ®P^)(f]  obtained  by  ueing  (0)  and 

(9)  ie  termed  the  bi linearly  blended  interpolant  of  F,  or  the  trensfinito  bi- 
*• 

linear  Interpolent  of  F.  Explicitly,  it  la  given  by 


U(e.t)  - (l-e)F(0,t)  ♦ sFtl.t)  ♦ tl-t)F(*,0)  ♦ tF(s,l) 

- (l-eMl-t)P(O.O)  - (1-e)  tF(0,l)  (10) 

- i(l-t)f(l,0)  - etra.l)  . 

*»  «♦ 

This  function  has  the  properties  that  U • P on  the  perimeter  of  the  unit  square 
[o.ljxto.l].  This  was  first  demonstrated  by  S.A.  Coons  in  (2). 

Figure  2 illustratee  the  mappings  induced  by  the  projectors  (6) -(8)  on  e 


•The  reader  should  verify  that  both  the  operators  P and  P are,  ir.  fact,  pro- 
jectors, i.e.,  they  are  linear  and  ideapotent. 
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region  R with  blending  functions  given  by  (9).  It  is  readily  seen  that  Pg[F] 

and  P [F]  each  match  only  two  opposing  boundary  segments,  P P [F]  matches  only 
t ^ 8 * 

the  comers  of  R,  but  (P  gP  )[F]  does,  in  fact,  interpolate  the  complete  peri- 

s ■* 

meter  of  R.  The  mappings  P ( F]  and  P [F]  are  sometimes  referred  to  as  "linear 

3 t 

lofting",  in  analogy  with  the  drafting  procedure  known  as  lofting.  P P [F], 

8 t 

which  is  linearly  ruled  in  both  directions,  is  termed  bilinear,  and  (P  « P )[F] 

s t 

is  properly  termed  bilinearly  blended.  The  three  projectors  Pg,  P and  Pg  e P 

+ 

are  all  transflnlte  projectors  since  they  interpolate  F at  a nondenumerable  num- 
ber of  points. 

Other  examples  of  transfinite  mappings  obtained  using  equation  (10)  are  shown 
in  Figures  3,4  and  5. 

We  can  generalize  the  above  notions  in  two  ways:  first,  we  may  consider  map- 

2 n 

pings  of  R ■*  E for  general  nj  and  secondly,  we  may  interpolate  F not  only  on 
the  boundary  of  the  region  R = { ,x^,  • . • ,x^)  - F(s,t):  0 s,  t <_  l} , but 

also  along  other  'flow  lines'  or  constant  generalized  co-ordinate  lines.  To 

< - 1,  and 


< 8,  * 1 and  0 ■ t.  < t,  < 
M 0 1 


thi3  end,  let  0 < s < sj  < 

M to 

let  {^(sO^q  and{^(t)}^oQ  be  functions  satisfying 

W “ 6ik'  VV  “ 5j i 

0<i,k<M,  0<j,  1 < H.  Now  define  the  projections 


(11) 


PS[FJ 


v 

1 $,(s)F(s  ,t) 

i-0 

N 

Pt[F]  5 I 

j«*0  3 


(12) 


The  product  projection 
H N 

p p (?)  s I I ♦1(8>*.tt)ru..t<>  (U) 

* * i-0  j»o  3 3 

interpolates  to  F on  the  finite  point  set  { ^ile  the  Boolean 
or  C rans  finite  interpoi ant 


SUB 


(P,  • ?t)(r)  s P(F)  ♦ Ptl?)  - P#PtlF] 


(14) 


interpolates  to  F along  the  H*N»2  lines  s » s^,  0<  1<  K and  t « t^ . 0 j N. 
That  is,  if  U(e.t)  ; IP  • Pt>[F)  then 


0(s,t^)  - K(s.t^) i 0 S j S S 
uts^.t)  - F(s.  ,t)  ,0SUK 


(15) 


If  H ■ tt  • 1 and  s^  * « 0,  * t^  <•  1,  (14)  reduces  to  the  transfinite 
bilinear  Interpolate  (10).  If  K * N • 2 and  »^  ” • 0,  e^  » - 1/2,  e^  " 
*2  «*  1,  then  using  the  blending  functions 
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$0(s)  - 2 (s-1/2)  (s-1), 

♦0(t>  - 2 (t— 1/2) (t-1) 

♦1(e)  - 4s(l-e) , 

^(t)  - 4t(l-t) 

(16) 

$2(a)  - 2s (s-1/2), 

4>2<t)  - 2t(t-l/2) 

in  (14)  yields  ths  blquedrstlcelly  blended  Intexpolent  of  F along  the  six  lines 
t - 0,  1/2,  1 and  t - 0,  1/2,  1. 

Thera  are  caaaa  in  which  bilinear  tranafinite  napping  techniquee  will  not 
produce  a aatiafactory  curvilinear  grid.  Thia  typically  happens  on  regions  R 
which  are  so  grossly  distorted  that  there  is  either  too  greet  a variation  in  the 
size  of  grid  elements  or  protions  of  the  generalized  co-ordinate  curves  actually 
aap  outside  the  region  ("overspill”,  cf.  [5]).  An  example  is  shown  in  Figure  6. 
Here,  a bilinear  transfinite  mapping  was  executed  with  the  result  that  some 
constant  generalized  co-ordinate  lines  overspilled  the  region. 

There  are  basically  three  ways  to  deal  with  such  difficulties.  First,  one 
may  decompose  the  overly  c cep lax  region  into  two  or  more  simpler  subregions 
and  aap  each  cf  these  separately.  Although  this  approach  generally  works  well, 
there  may  be  problems  at  the  interfaces  between  subregions,  since  the  general- 
ised co-ordinate  lines  will  have  slope  discontinuities  there.  One  way  of 
handling  this  difficulty  is  to  employ  higher  degree  blending  functions,  e.g., 
cubic  Hermit*  blending  functions. 

A second  way  of  attempting  to  achieve  a satisfactory  transformation  is  to 
reparametrize  the  boundary  segments  of  R by,  for  example,  introducing  monotonia 
transformations  of  the  independent  variables  e and/or  t. 

Another  way  of  dealing  with  complex  regions  is  to  introduce  auxiliary  con- 
straints into  the  transformation  problem.  Since  the  paramount  objective  is  to 
obtain  a one-to-one  (invertible)  mapping  of  S onto  R,  the  analyst  is  perfectly 
free  to  enforce  whatever  additional  cons  -ai-ta  he  feaie  will  guarantee  the  in- 


vertibility  of  the  mapping,  in  our  expe»*«  :>■  , j have  generally  found  it  ads- 
quate  to  specify,  ee  an  auxiliary  constraint,  the  image  (i  s.,  the  mapped  posi- 
tion) of  e single  Interior  point  of  S.  That  is,  v»  identify  where  inside  R we 
desire  to  sap  e selected  point  in  the  interior  of  $.  For  simplicity,  suppose 
the  point  in  S whose  image  position  we  want  to  constrain  is  the  mid-point  of 
the  square,  a « t - 1/2.  We  went  to  force  this  point  to  nap  into  the  point  in  R 
having  co-ordinate*  (a.B) • Let  U be  the  bilinearly  blended  function  of  (10). 
Then,  the  following  transformation  maps  2S  onto  3R  and  maps  (1/2, 1/2)  onto  the 
point  (a  ,0)  > 

$<s,t)  - 5(e,t)  ♦ 16e(l-s) t(l-t) 


- 11(1/2,1/2) 


(17) 


To  verify  this, note  that  along  the  perimeter  of  {0,l)s(0,l]  the  right-hand  side 


t77 


reduce*  to  just  t5 (s,t) , which  naps  35  exactly  onto  3ft,  as  desired.  Por 
(*,t)  ■ (1/2, 1/2),  the  right-hsnd  side  reduces  to  just  (a, 6) . More  generally, 
the  point  (s,t)  « (a,b)  can  be  napped  into  (0,0)  by.  the  foraulat 

[(»)  • 5<,'t>]  • 1181 

As  an  exanple,  Figure  6 shows  a region  R for  which  the  bilinear ly  blended 

transformation  (10)  does  not  give  an  invertible  mapping  of  5 onto  R.  It  is 
intuitively  clear  that  the  inage  of  the  point  (s,t)  ■ (1/2, 1/2)  has  mapped  too 
far  to  the  right.  Therefore,  we  enforce  the  auxiliary  constraint  that  the  point 
(1/2, 1/2)  in  5 should  map  onto  the  point  (.494,. 119)  in  R.  Figure  7 illustrates 
the  result  of  the  transformation  obtained  using  (17) . 

In  [5],  Gordon  and  Hall  propose  the  use  of  more  general  auxiliary  constraints. 
For  example,  instead  of  just  a single  point,  they  consider  mappings  for  which 
lines  of  constant  s or  t are  forced  to  map  onto  preselected  generalised 
co-ordinate  curves  in  the  domain  R.  Such  curves  may  arise  naturally  cs  inter- 
faces between  subregions  of  R,  or  they  may  be  determined  by  the  analyst  on  geo- 
metric grounds.  The  transformation  formulas  appropriate  to  these  constrained 
maps  are  given  by  aquation  (14) . 

If  the  region  R is  basically  triangular,  in  contrast  to  quadrilateral  trans- 
finite  interpolation  techniques  over  triangles  may  be  more  appropriate.  The 
theory  for  such  "trilinearly  blended*  methods  was  developed  in  (1).  The  details 
of  these  techniques  aa  applied  to  grid  generation  may  be  found  in  [6]  and  (7). 

As  a practical  matter,  the  curvee  bounding  R may  not  be  easily  repreaentad 
as  closad-foxm  mathematical  expressions.  Nevertheless,  the  above  results  still 
apply  if  the  boundary  curves  are  represented  es  discrete  point  sets,  i.e., 
piecewise  linear  curves.  For  e fuller  discussion  of  “discredited  t nuts  finite 
mappings"  tee  (6)  and  (7].  Surfaces  in  Euclidean  3 -space  are  handled  in  pre- 
cisely the  mm  way  as  2-0  regions.  All  that  nead  be  dona  la  to  add  the  third 
co-ordinate  functions  S(e)  end  Z(t)  to  the  x end  y components.  A discussion  of 
surface  decomposition  techniques  is  given  in  {7]. 

3.  THAEE-DIHEKSIOKAL  SOLID  S7WCTUSE5 

The  purpose  of  this  section  la  to  outline  the  extensions  to  1 -dimensions  of 
the  bivariate  tranefinlte  mapping  techniques  discussed  above.  To  begin,  we 
consider  the  following  three  projectors t 
PjF)  - (l-*)F(0,t,u)  ♦ sFd.t.u) 

r)  - (l-t)F(s.o.u)  ♦ tF(s,l,u) 

• (l-u)F(*,t,Q)  ♦ uF(s,t,l)  . 


(19) 
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In  these  expressions,  the  "primitive  function"  t is  4 vector-valued  function 
of  the  three  independent  parameters  s,t  and  u.  As  a,t  and  u range  over  the  unit 

cube,  F maps  out  the  solid  volume  R.  Depending  upon  the  co-ordinate  system  used, 

*♦> 

the  three  components  of  F may  correspond  to  Cartesian,  spherical,  cylindrical, 
toroidal,  etc.  co-ordinates.  As  a practical  matter,  the  co-ordinate  system 
employed  should  be  that  which  most  appropriately  fits  the  geometry  and  topology 
of  the  problem  domain. 

Quite  clearly,  the  three  projectors  in  (19)  correspond,  respectively,  to 
linear  blending  (lofting)  in  s,t  and  u.  Now,  however,  the  geometric  entities 
on  the  right-hand  side  of  the  expressions  are  not  curves,  as  in  (6),  but  rather 
surfaces.  For  example,  as  t and  u range  over  the  parameter  domain  [0,l]x(0,l], 
the  vector-valued  function  F(0,t,u)  traces  out  a surface  in  Euclidean  3-space 
corresponding  to  one  of  the  six  faces  of  the  solid  volume  R. 

It  is  useful  to  consider  the  pairwise  products  of  the  above  three  projectors. 
For  instance,  the  product  of  the  first  and  the  second  is 

P (P  [F])  “ (1-s)  (l-t)F(0,0,u)  ♦ (l-s'tE(0,l,u) 

3 t -*  -»  (JO) 

♦ s (l-t)F (1,0 ,u)  ♦ 2tF(l,l,u) 

The  right-hand  side  of  this  expression  contains  four  expressions  which  refer  to 

the  edges  of  the  object  under  consideration.  By  evaluating  the  expression 

along  the  four  edges  (s,t)  - (0,0) , (0,1) , (1,0)  and  (1,1),  it  can  be  verified 

that  the  tri variate  function  P P fF}  does,  in  fact,  match  F along  these  edges. 

s t .*  ♦ 

(It  m.  , easily  demonstrated  that  P P [F)  • P P [P],  i.e.,  the  projectors 

® t t 3 

commute,  just  as  in  the  bivariate  case.) 

We  have  seen  that  the  projectors  Pg,Pt  and  Py  each  interpolate  the  two  oppos- 
ing faces  of  the  solid  volume  described  by  F(e,t,u),  and  that  products  of  pairs 
of  these  projectors  interpolate  to  tins  edges  of  the  solid.  Xf  we  take  the  pro- 
duct of  all  three  of  the  projectors  in  (19),  we  obtain  the  expression! 

P P P (?)  - (1-s)  (1-t)  (l-u)P  (0,0,0)  ♦ (1-s)  (l-t)uf(0,0.1) 
a t u 

♦ (1-s) t(l-u)F (0,1,0)  ♦ (l-#)tu?(0,l,l) 

* * 

♦ e(l-t)  (l-u)f  (i.0.0)  ♦ s(--t)uP(l»0,l) 

♦ stu-siru.i.o)  ♦ stum. i,i>. 

The  right-hand  side  of  this  expression  contains  values  of  F which  refer  to  the 
eight  corners  of  the  region  S.  It  does,  in  fset,  interpolate  to  these  eight 
ooroera.  Since  (21)  is  linear  ir.  each  of  the  three  parameters  s.t  and  u,  it  is 
termed  a trillnear  Ineerpolan: . The  graph  of  the  trilinear  interpolant  ie  a six- 
sided  polyhedron  which  passes  through  the  vertices  of  R»  it  is  simply  the  3-0 
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generalization  of  a quadrilateral  in  Euclidean  2-space. 

In  the  2-dimensional  case,  one  starts  with  the  two  projectors  P^  and  P of 
(6)  and,  by  combination,  generates  a total  of  four;  namely,  P , P , P P and 

S t 9 t 

P • P^.  In  three  dimensions,  the  situation  is  much  more  complex  and  the  variety 

of  possible  projectors  much  richer.  In  [3]  and  (4],  it  is  shown  that  under  the 

two  binarary  operations  of  operator  multiplication  and  Boolean  (•)  addition,  any 

collection  of  commutative  projectors  forms  a distributive  lattice.  Space  does 

not  permit  going  into  the  details  of  this  theory,  but  we  can  illustrate  some  of 

the  results  in  the  trivariate  case.  Without  proof,  we  state  that  there  are  21 

distinct  projectors  which  can  be  formed  by  multiplication  and  ® addition  of  the 

three  projectors  P , P and  P . In  addition  to  those  displayed  above,  the 
s c u 

following  are  examples: 


Ps* 


P ♦ P - P P„ 

S t St 


P 9 PP  * P + PP  - P PJ> 
s tu  a tu  stu 


P P • P 

St  u 

p • p 
a t 


P P ♦ P - P P P 
st  u stu 


PP 
t u 


p ♦ p - p P 

S t St 


PP  9 ° V 9 P ? 

J t t u US 


PP  * PP  * P P -2PPP 
st  tu  us  stu 


Ps*Pt 


p ♦ p ♦ p - PP  - PP 
s t u st  t u 


P P ♦ P P P 
us  stu 


(22) 


One  should  note  that,  because  of  the  ideapotency  and  linearity  of  the  projectors, 

much  cancellation  occurs.  For  instance,  one  has  P ♦ ? F * P ♦PP  -PPP 

• st  s at  b s t 

• P , which  means  that  nothing  is  achieved  by  taking  the  • sum  of  P and  P P . 
this  is  because  the  interpolation  properties  of  P#pt  are  a subset  of  thoso  of 
the  projector  P . 

an  aspect  of  the  theory  developed  in  {3}  and  {4]  ie  that  there  is  an  isomt- 
pMsm  between  the  distributive  lattice  of  projectors  and  the  associated  distri- 
butive lattice  of  their  precision  sets.  In  other  words,  if  wo  know  the  expres- 
sion for  a certain  projector,  then  we  can  determine  the  point  set  on  which  it 
interpolates  by  replacing  operator  multiplication  by  sat  intersection  and  • 
addition  by  set  union.  For  example,  the  precision  set  (set  of  points  on  which 
it  interpolates)  of  P#  consists  of  the  two  faces  of  ft  defined  by  s • 0 and 
a » It  and  similarly  for  P and  Py.  thus,  it  follows  from  the  isomorphism  that 
the  precision  sate  of  the  projectors  in  (22)  are,  respectively,  given  by  the 
following  expressions  in  which  s . 5.  and  5 are  the  precision  sets  of  P , P, 

ft  v U 9 s 

aid  p i 
U 
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S u S_ 

■ t 

S u (S  n S ) 
s t u 

(S  o S )uS 
8 t U 


S u S.  u (S  n S ) 

• t t u 

(S  n S ) u (S  n S } u (S  r>  S ) 
st  t u us 

S u S u S 

8 t U 


(23) 


If  one  thinks  of  the  precision  sets,  it  is  obvious  thet  the  weakest  (‘algebrai- 
cally minimal'')  of  all  projectors  is  the  triple  product  P^P^P^  end  the  “algebrai- 
cally maximal"  projector  is  the  Boolean  sum  of  ell  threat  P#  • Pfc  • P . XI 1 
other  possible  projectors  are  algebraically  in-between  these  two. 

Which  of  this  airiad  of  3-D  projectors  one  uses  in  practice  ie  a latter  of 
what  data  is  given;  or,  wore  precisely,  where  the  data  is  given.  In  other  words, 
one  is  given  the  precision  set  and  mist  use  the  isomorphism  ‘backward"  to  infer 
the  appropriate  projector  (i.e.,  interpolation  formula).  For  instance,  if  one 
does,  in  fact,  know  the  exact  shapes  of  all  six  of  the  bounding  surfaces  of  R, 
then  the  appropriate  mapping  equation  is  the  full-blown  expression 
(P#  e e Pu)(F],  which  explicitly  Involves  all  faces,  ell  edgas  and  all 
corners.  Xt  the  other  extreme,  one  nay  only  know  the  co-ordinates  of  the  eight 
corners  of  K.  In  this  case,  one  would  use  the  transformation  equation  (31). 

In  practice,  the  situation  ia  usually  somewhere  in-between,  i.e.,  the  date  ie 

aeldom  ee  casplete  ee  required  by  the  maximal  projector  P • P • P or  es  meant 

« a t u 

aa  only  the  eight  vertices  needed  in  P^P^P^lF] . The  examples  given  beltw  assume 
that  the  function  r ie  known  (i.e.,  data  ie  given)  on  the  13  edges  of  R.  In 
this  ceea,  the  revelent  transformation  equation  or  mapping  formula  is> 


(P  P • P P IP?  )[f) 
at  t u u a 

- (1-a)  (l-t)F(0.0,u)  ♦ (l-e)trtO.l.u) 

♦ a(i-t)?Cl,6,o*  * st£(l,l,u) 

♦ (1-t) (l-u)?(s.0,0)e(l-t)u?(a,0.1) 

♦ t(l~u)?(s,l,0)  * tu^(a.l.l) 

♦ U-aHl-u)P(O.t.O)  ♦ a-s) u£(0,t,l>  (34) 

♦ e(l-u)?U.t,0)  » suFU.t.l) 

- 3(  (1-a)  (i-t)  (1-u)  ?(0,0,0)  * (l-s)(l-t)uirC0,0,i) 

♦ U-«)t(l-u)?(0,1.0)  * (i-s)t«r(0,i,l) 

♦ a(l-t)  U-u)?()  ,0,0)  * eU-t)u?!l,0.1> 

♦ st(l-o>?(1.1.0)  * etu?U,  1,1)1. 


The  curvilinear  cc-ordlnate  system  generated  by  the  following 
paranetrlzatlar. 


F(0,t)  = 
F(s,0)  = 

FvS#l)  = 


0.5t  + 2 
-2 

'1,25  - .75  COS  (rs) 

-2  + ,75  sin  tfTS)_ 


-2.5t  + ,5 

“2 

, F(l,t)  = 

<-6,5  + 15, 5s, 
2,5, 

V 

-2  + 6s, 

- 0,0, 

4 - 6s, 


0,0  < s < .33 
,33  < s < .66 

,66  £ S £ 1 

0,0  £ s < .33 
.33  £ s < ,66 

,66  £ S £ 1 
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Figure  9 

TTk  transflnlte  m®  of  a regicn  with  a crack  via  the  following 
paranetrlzatlcn : 


F<0,t)  = 


*1  - 3t 
-3 


F(l,t)  “ 


*1.05  + l,45t 
-3 


F(S,Q)  = 


F<$,1>  = 


1 + .05s 

-3  ♦ 3s,  0.0^5  10.5 
-3$,  0.5  i s 1 1 , 


-2, 

O.Q  < S < .33 

-12.7  4-111  ,<is  -27.9s1 2, 

.53  < s < .66 

2.5, 

.66  < S < 1 

-3  ♦ 9s 

0.0  < s < .33 

0.0 

.33  < s<  .66 

6 - 9s 

.66  i s i 1 

L. 


Figure  5 

curvilinear  co-ordinate  systan 
TMe  bandary  segments  are: 


F(0,t) 


F(l,t)  = 


W. 


F($,0) 


**  - 11.5s  + 7,5s2,  0.0  < s < ,33 
’ q-»s  ^ 36s2.  ,33<  s<  .66 
[-3.5s  + 7.5S2,  *66  < s < I 

f-1 , 125s  ♦ l , 125s2,  0,0  <,  s < ,33 

I *75  + ],5s  ,33  < s < ,66 

1.I2&  - 1.125s2  <66  < S£  i 

*»  - 20s  ♦ 20s2  * 

L'2  + 4s  J , 


F(s,l) 
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Figure  6 

Blllnecr  blending  here  yields  a non-mlvalent  man,  The 
Dararetrlzatlcn  of  tie  boundary  Is: 


0 

'l  + 2t  - 212 

R0,t)  » 

F(U>  3 

_t. 

X 

's' 

-*• 

r 

s 

F(s*0)  = 

, F(b>  1 ) «* 

\ - 3s  ♦ 3s2. 
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Figure  7 

A curvilinear  co-ordinate  system  with  no  "oversell  1"  is  achieved 
via  (17).  Here,  the  point  (1/2,172)  in  the  s,t-olcne  isnroped 
onto  the  point  (.494,  .119)  in  the  x,y-olcre. 
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’ w v 


ZnTlZ  V'°  a8)  * W*B  the 

^ V®  ln  «*  S't-olone  mto  the  ootnt 
(•JS-.33W  !n  the  x,H)l(ne.  The  bomry  segrents  ore: 


™'()”  [JJ-  F«-t)  . 


-3+6  cos 
-3+6  sin 


F(s,0) 


6s  - 3 


F(s,i) 


|~3  ♦ yj~^5~~  (6s  - <J.5)2  , 
1 ^ * 

|-3  * /.25  - (6s  - 0.5)2  , 
[,6s-  3 


(trt/2) 

.67  i s < .83 
otherwise 

M 

• 

•67  £ s i .83 
otherwise 
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